Aharonov-Casher Effect for Plasmons in a Ring of Josephson Junctions 
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Phase slips in a one-dimensional closed array of Josephson junctions hybridize the persistent 
current states (PCS) and plasmon branches of excitations. The interference between phase slips 
passing through different junctions of the array makes the hybridization sensitive to the charges 
of the superconducting islands comprising the array. This in turn results in the Aharonov-Casher 
effect for plasmons, which in absence of phase slips are insensitive to island charges. 



For over two decades, arrays of small superconduct- 
ing islands connected by Josephson junctions have been 
a popular platform for experimental studies of quan- 
tum fluctuations in superconductors |l|. This popu- 
larity has originated in part from the tunability of the 
islands' charging energy and the junctions' Josephson en- 
ergy, which has enabled the design of arrays with a de- 
sired level of quantum fluctuations of the phases of the 
order parameter. 

The most fundamental manifestation of the quantum 
fluctuations of phase is a quantum phase slip. The pro- 
liferation of phase slips in long nanowires or Joseph- 
son junction arrays is predicted to lead to a quantum 
phase transition between superconducting and insulating 
states |2j, y] ■ In addition, coherent phase slips interfere 
with each other, giving rise to the Aharonov-Casher ef- 
fect [4[ . An unambiguous observation of these predictions 
in dc measurements has proven to be difficult |l|, la, [fjj ■ 
However, an evidence of the Aharonov-Casher effect in 
long arrays has recently been seen in spectroscopic mea- 
surements [3] performed on a ring of junctions pierced by 
a magnetic flux. The measurements [3, |8| focused on the 
avoided crossing of the two lowest-energy states. These 
two states carry counter-propagating persistent currents 
and become degenerate if the magnetic flux equals to 
a half-integer of flux quantum; phase slips remove the 
degeneracy. The Aharonov-Casher effect reveals itself 
through the modulation of the avoided-crossing gap by 
charges in the superconducting islands. To interpret the 
experiment of Ref. [3] , it was sufficient to account for the 
hybridization of two otherwise degenerate lowest-energy 
many-body state s M - I ffi - |l Ol| . However, measurements of 
the same type [3, la, |ll|, LL2J routinely display a rich, flux- 
dependent level structure at higher energies. This struc- 
ture is poorly understood and, to our knowledge, the 
effect of phase slips on the spectrum of excitations has 
not been studied theoretically. 

In this work, we investigate the effect of quantum phase 
slips on the spectrum of plasmons in a one-dimensional 
array of Josephson junctions forming a ring. In a dense 
plasmon spectrum, the structure of the multiple avoided 
crossings and their sensitivity to the charges (Aharonov- 



Casher effect) turns out to be quite different from the 
single avoided crossing |9|, |lfj] occurring in the ground 
state, see Figs. [2] and [3] 

We consider a ring of TV ^ 1 identical junctions (see 
Fig. [I]) of Josephson energy Ej. The phase (pi of the 
order parameter may be taken [l| uniform within the 
grain i. Each grain is capacitively coupled to a volt- 
age source V* with gate capacitance Co; this models in- 
evitable [3] quasi-static charges induced by the environ- 
ment. Two neighboring grains are connected through a 
Josephson junction shunted by a capacitance C. The ring 
is threaded by a magnetic flux $ that induces a persis- 
tent current |lfj] . Since the self-inductance of the ring is 
negligible, 3> is close to the flux of the externally applied 
magnetic field. 
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FIG. 1: Segment of a ring of Josephson junctions. Grey circles 
represent superconducting islands, which are connected by 
ideal Josephson junctions (crosses). 

At Vi = 0, the Hamiltonian of the ring is 



N 



N 



n = 2e 2 Y, Q l {C- 1 ) lJ Q J + Ej J2 [1 - cos (7*)] , (1) 

where the first and second terms are the kinetic and 
potential energies describing the energy cost of elec- 
trostatic charging and persistent currents, respectively. 
Cij = 5i,j(2C+Co) — (Sij+i +d~i : j-i) C is the capacitance 
matrix, and ji = <Pi+\ — (p% + S/N are the gauge- invariant 
phase differences across the junctions, with 5 = 27r$/"I>o. 
The operator Qi corresponds to the charge on island i 
in units of — 2e and is canonically conjugate to (pi. The 



topology of the ring enforces the constraint [lc 



N 

8 = 1 



(5 = (mod 2tt). 



(2) 



Equation ([1} is supplemented by the boundary condition 
tp((pi,...,(fi +2TT,...,if N ) = il>((pi,...,(pi,...,(pN), for 
all i € {1, ... , iV}. Because only phase differences enter 
Eq. (|Tj). the problem is invariant under a simultaneous 
shift of all (pi) the associated conserved quantity is the 
total charge Q = Yl < Q« G Z. Hereafter we restrict the 
discussion to a fixed value of Q. 

When V~i ^ 0, induced charges g, = C Vi/2e break the 
rotational symmetry of the array and result in a shift of 
the charge operators in Eq. ([1}: Qi —t Qi — qi- A unitary 
transformation eliminates the ^-dependence from 7-L and 
makes the wave functions independent of the total phase 
4> = ^j (fij. The outcome of this transformation is a 
Hamiltonian given by Eq. JT]) (modulo constant terms), 
and a boundary condition [9] 



ip(<Pi ...,<Pj+2n,..., <p N ) 



,^,...,^jv). (3) 



The main objective of this work is to evaluate the 
energy spectrum for Eq. ([TJ. We are interested in the 
regime of Ej dominating over the charging energy. To 
zcroth order in the charging energy, the phase variables 
become classical. The configurations minimizing the po- 
tential energy in Eq. (JXJ) under the constraint of Eq. @ 
are characterized by an integer m, 



ji mod 2ir = (27TTO + 6)/N . 



(4) 



for all i E {1, . . . , N}. These states, which host an elec- 
trical current whenever $ ^ (mod $o)i will be referred 
to as "persistent current states" (PCS) \m). The specific 
value of m leading to ji mod 2n £ [—n/N, n/N] describes 
the classical ground state, while other values of m denote 
the local minima of the potential and correspond to ex- 
cited (metastable) states with higher currents. 

A small but nonzero charging energy results in quan- 
tum fluctuations of % around the (local) minima. For 
low-current states (7, mod 2n ~ 1/JV <C 1), small fluctu- 
ations are described by quadratic expansion [14] . 



N-l 



Uo =J2 E m\m)(m\ + J2 m i U 
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(5) 



of Eq. ©. Here E m = (2nm + 5) 2 Ej/(2N), and a\ is 
the creation operator for a harmonic mode ( "plasmon" ) 
of frequency $7;. Note that only N — 1 modes appear 
in Ho [15]. Herein we restrict the discussion to zero- 
and one-plasmon states only, for which anharinonic ef- 
fects [111 are small. In addition, we will be interested 



in low-current-carrying states for which the plasmon fre- 
quencies are essentially independent of the magnetic flux: 
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where tu p and Eq = e 2 /2C are the single-junction plasma 
frequency and charging energy, respectively. To sim- 
plify notation, we have labeled the plasmon modes by 
I € {1, ■ • ■ ,N — 1} instead symmetric about I = wave 
numbers [16J. At Co ^ 0, one-plasmon states are doubly 
degenerate (except for I = N — I) because £li = fijv-z- 
At Co — 0, the degeneracy increases to N — 1 because 
£li = Lu p for all I (dispersionless plasmons) . 

The harmonic model Ho provides a good approxima- 
tion for the energy spectrum of the ring, except in the 
vicinity of flux values for which two neighboring PCS 
(\m) and \m — 1)) become degenerate. These degenera- 
cies, which emerge at $ = $o/2 (mod $o)> can be lifted 
by tunneling. Tunneling between \m) and \m — 1) in- 
volves a change of ~ 2ir in the phase difference across 
a junction, and is known as a quantum phase slip. We 
neglect tunneling processes involving multiple simultane- 
ous phase slips, because they are highly unlikely in our 
regime of interest (Ej >• Ec)- 

The influence of phase-slips in the energy spectrum of 
"Ho can be incorporated perturbatively by adding 5T-i^ 
and 5%^ in the zero- and one-plasmon subspaces, re- 
spectively, where 

sn {0) = J2 (^ (0) \ m - 1) H + c - c -) 

m 

N-l 

<^ (1) =E E (v$\rn-l,l)(l',m\ + c.c.). (7) 
m ;,i'=i 

In Eq. (J7J), i/°) is the \m) — > \m— 1) tunneling amplitude 
in absence of plasmons, whereas v\ J is the \m) — > \m— 1) 
tunneling amplitude in presence of one plasmon that gets 
scattered from mode I' to /. For $ ~ $o/2 (mod $ ) 
and Co — 0, the splitting between the two lowest-current 
levels is 



8E^=[{E m -E m _ x f + A\u^ 
8E^ = [{E m -E m _ 1 f + A\n n \] l/2 



1/2 



(8) 



in the zero- and one-plasmon subspaces, respectively. In 
Eq. ©, jU„ is the n-th eigenvalue of the matrix i/Wtj/W 
(n <G { 1 , . . . , A — 1}). The flux- dependence of the splitting- 
comes mainly from £"„,.($), while v^°> and /i„ can be 
evaluated at $ = $o/2. 

Next, we proceed with the microscopic derivation of 
the tunneling amplitudes. While the expression for ia°) 
is well-known, this Letter develops the first theory for 
j/ 1 ) p or convenience we work with the non-gauge- 
invariant phase differences 6i = tpi+i — fi- Without 



loss of generality, we consider the two lowest-current 
states at <& = <&o/2 (\m = 0) and \m = -1))- Each of 
them is represented in #,j-space by a manifold of points 
that minimize the potential energy [Eq. (|3|)]. Represent- 
ing \m = 0) by {9i} = 0, its N neighboring minima 
{9i} = -2ir/N + 6i, n 2Tr (where n E {1,...,N} all repre- 
sent the \m = —1) PCS. Starting from {6^} = 0, each of 
the \m = —1) minima can be reached through an instan- 
ton on a least-action path. Along each of these paths a 
different junction (labeled by n) undergoes a phase slip. 
The total \m = 0) — > \rn = —1) tunneling amplitudes are 
the sums over the contributions of N individual paths 



(denoted by v n and ^„.j ;/) 
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and 
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(9) 



Plasmons enhance the tunneling amplitudes via the os- 
cillation energy available to assist a phase slip. Different 
junction variables oscillate with different amplitudes, de- 
pending on which of the TV — 1 plasmon modes is excited. 
That is why contributions of different junctions are differ- 
ent even as the instanton paths end in equivalent points. 

The calculation of v\ J requires specifying the plasmon 



wave functions. For \m — 0) and Co 
chosen as 



0, they can be 



N-l 



in =KtJ2 °ne l 



,-D*(j: N n r = \e n ,) 2 



N-l 

n< 
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where I € {1, . . . ,N - 1}, D 2 = hu p /(16E c ), and K x oc 
[1 + ^i,n-i{N x / 2 — 1)] is the normalization factor. For 
the N neighbouring minima representing \m = —1), we 
can pick exactly the same wave functions up to a shift in 
coordinate space, and denote them fy . Due to Eq. ([3]), 
■0; and tpl are related by a phase factor. 

The next step is to find the contribution of a single 
junction n, namely v n \ v . Three difficulties become ap- 
parent. First, a phase slip occurring in junction n in- 
volves not only 9 n but also requires a shift of —2ir/N in all 
9i^ n - Second, a harmonic approximation of the potential 
energy of junction n is no longer appropriate because j n 
changes by almost 2ir. Third, having N — 1 modes for N 
junctions means that plasmons are delocalized along the 
ring, which in turn requires a prescription to determine 
how a plasmon assists the phase slip in junction n. We re- 
solve the first two difficulties by (a) changing variables by 
a linear transformation {9i} i-> {xi} [l7| such that only 
x\ changes during a phase slip, and (b) keeping the co- 
sine term for the n-th junction while making a harmonic 
approximation for the rest. The minimum {9{\ = corre- 
sponds to {xi} = and {9i} = —2ir / N + 5i^ n 2ir becomes 
{x^ — 5i : i2ir(N — 1)/N. Importantly, the proper choice 



of {x^ simplifies Eq. (fTJ) into H e & + Hharm, where 

1 



Wett — Ej 

AE C N - 1 

+ ~h 2 n~ 



l -^ {xi+ w)-^ xi+ m^~f i2 



Px] 



is an effective Hamiltonian for junction n, and 



N-l 

%harm = / J fu^p ( — 

k=1 V " 



b{h 



Ej N-1 

2N N ' 



(10) 



(11) 



describes N — 2 harmonic modes that are decoupled from 
x\. The operator p Xl in Eq. (|10[) is conjugate to x\, and 
the linear and quadratic terms in x\ arise from Eq. (fTJ) 
and Eq. ©. 

A harmonic approximation of "H c ff + "Hharm about a 
minimum results in a set of modes {6fc}, which are gen- 
erally different from {fife} of Eq. |5]). In terms of these 
new modes, only 6}, associated with variable X\, creates 
an excitation affecting tunneling in junction n: in absence 
of this plasmon the amplitude has some absolute value 
z/o, and if there is a plasmon the amplitude has a different 
absolute value v\ [181 ]. Thus, the initially multidimen- 
sional tunneling problem is now reduced to a textbook 
one-dimensional problem. Next, we expand ipi> and ip} n 
in terms of the two sets of plasmon wave functions that 
follow from the harmonic approximation of H e e + Hhaim 
about x\ = and X\ — 2tt(N — 1)/N, respectively. This 
resolves the third difficulty mentioned above. From the 
expansion, we arrive at 



,(!) 



2™(l-,' 



V n<N;l,l> = Zn v Q$l,V + a l a l' e% N ~ 1 z n{v\ - Vq) 

v n]i,v ~ S W N + Sln-i(vi - V )] , 
where a t = [N 1 / 2 + Si, N -i(l - N^ 2 )]/(N - 1), and 



(12) 



z n = exp 



n / , N 



Qk 



fe=i 



* + iV 



(13) 



is a phase factor (defined up to a global phase) that arises 
from Eq. (J3J. In Eq. (TT2"]) . the second term in the right 
hand side describes the enhancement of tunneling due to 
plasmon oscillations having a nonzero component along 
the direction of the phase slip. A simplified version of the 
procedure described above can be applied to recover [9( 
the ground-state tunneling amplitude, r/„ = z n Va- 

With no induced charges, the tunneling amplitudes are 
sensitive to the total charge Q. If Q/N £ Z, we obtain 
[see Eqs. ([8]l, (|9]l, (pj). and (p])] i/°) = Nv and 



vft, =8i, v {Nu Q + e), 



{v\ - vo) 



N 



N-l' 



(14) 



i.e., n n = (NvQ + e) 2 . Instead, if Q/N ^ Z it follows that 



(0) = |3|, /ii = and /U„>i = e 2 . Thus, destructive 
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FIG. 2: Conditional distribution of splittings of the one- 
plasmon level for given |i^°'|, according to Eq. (|16[1 . Shown 
is the case e < Iz^ '!. 



interference caused by Q/N ^ Z protects a 2-fold degen- 
erate crossing point. The one-plasmon energy spectrum 
remains highly degenerate at any Q/N. 

In reality, there are induced charges and there is no 
control over them [7|. Below we regard qi as indepen- 
dent and identically distributed (i.i.d.) random variables, 
and model them by Gaussian probability density func- 
tions (PDF) with mean zero and standard deviation a. 
Consequently, z n and the tunneling amplitudes become 
random variables. For broad distributions (<j > 1), the 
variables z n are approximately i.i.d. with uniformly dis- 
tributed phase, and the tunneling amplitudes become in- 
sensitive to the value of Q. 

At N ^$> 1, the PDF for the ground state tunneling 
amplitude follows from the central limit theorem, 



V{\v 



(0)1 



2|i/°>| 



; • '-n- 



"/(, 



2|i/°>| 
Nv 



(15) 



where I n (x) is the modified Bessel function of the first 
kind. From Eq. (fl"5j) . we obtain (to leading order in N) 
the expectation value (|i^|) — uo>/Ny/n/2. 

For the plasmon-assisted tunneling at N ^> 1 we 



find [17J Hn = \v 



(0)12. 



-T) n , where r\ n sa 2eRe(^°)z* 
Thus the conditional PDF of \^ n \ 1 ^ 2 for a given \u^\, 



V{\Hn 



V2||,/0)|\ 



D : 



|M» 



11/2 



n y/4e 2 \v(°)\ 2 - (\ t i n \ - \i>W\ 2 - e 2 y 



|I/(0)|2_ e - 

(16) 

is finite in the window \\u^\ — e| < l^nl 1 ^ 2 < k (0) l + e , 
see Fig. H When |i/°)| > e, the distribution of Eq. ([16]) 
is approximately valid for all eigenvalues /i n . When e > 
|z^°)|, there is one eigenvalue that does not conform to 
Eq. ([T5)l . This single eigenvalue is smaller than |i/°) — e|, 
and corresponds to a spatially delocalized mode [171, |20j . 
The preceding paragraph evidences that the 2(N — 1)- 
fold degeneracy in the one-plasmon subspace of Ho is 
completely removed by phase slips in presence of ran- 
dom charges (see Fig. \S§. This is a manifestation of 
the Aharonov-Casher effect for plasmons. Inverting the 
unitary transformation that approximately diagonalizes 
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FIG. 3: Schematic energy spectrum for a ring of JV Josephson 
junctions (Fig. [TJ with Josephson energy Ej, small ground 
capacitance Co, and a given realization of random gate- 
induced charges. The ring is pierced by a magnetic flux 
$. (a) Spectrum in absence of plasmon excitations, where 
E m = ty 2 Ej/(2N) + (N - l)liuj p /2 and uj p is the plasmon 
frequency, (b) Spectrum in presence of one plasmon, where 
E (i) = E (p) + ^ In thig figure |i/C0) | > £> The non _ 

degeneracy of the 2(N — 1) one-plasmon energy levels is due 
to the Aharonov-Casher effect. The dashed lines illustrate the 
energy spectrum in absence of phase slips. 



I/ ( 1 )tjy(i) (which turns out to diagonalize i/W as well), 
and applying it onto the plasmon states ipi, we obtain 
N — 2 modes (N — 1, if |f'°)| > e) that are almost per- 
fectly localized in space. Hence, in presence of randomly 
induced charges, phase slips lead to the localization of 
all-but-one (all, if |iA°)| > e) plasmon excitations in sin- 
gle junctions. Without phase slips, the plasmon modes 
would be insensitive to induced charges. 

Thus far we have assumed dispersionless plasmons, 
21|. For Co 7^ 0, the results derived above still 



C o ~0 

hold provided that the bandwith of the plasmon disper- 
sion (ojp — fli=i) is small compared to e. In the opposite 
case, N 2 Cq/C > 16tt 2 e/(Hu) p ), the effect of C on the 
plasmon energies can no longer be neglected. For large 
enough Co , the typical energy separation between neigh- 
boring plasmon modes exceeds e, and the hybridization 
due to phase slips is restricted to subspaces of four-fold 
degenerate modes (two PCS and two plasmon modes I 
and N - I). As long as N 2 C /C < 47r 2 , the tunnel- 
ing amplitudes between these four levels can be found by 
the method developed above, but with a set of plasmon 
modes diagonalizing Eq. ([5]) for Co 7^ 0. In this case too, 
tunneling in presence of random induced charges removes 
all degeneracies. Yet, unlike for Co = 0, plasmons remain 
delocalized, thus weakening the Aharonov-Casher effect. 
In summary, we have presented a theory for quantum 
phase slips in a closed Josephson junction array in a one- 
plasmon excited state. Phase slips hybridize different 
one-plasmon and persistent current states, and lead to 
a typically non-degenerate energy spectrum that is sen- 
sitive to gate-induced charges. This sensitivity is a sig- 
nature of the Aharonov-Casher effect. Future work will 
address higher excited plasmon states, the interconver- 
tion between current and plasmon excitations, and the 



effect of thermally excited quasiparticles. 
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This supplemental material has two parts. First, we 
provide technical details for the derivation of Eq. (10) and 
Eq. (11) in the main text. Second, we provide technical 
details that justify Eq. (16) in the main text, and identify 
a single eigenvalue that can depart from it. Note about 
references: the references cited below are listed at the 
end of the supplemental material. 
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In this section we provide explicit expressions for the 
coordinate transformations that result in a decoupled ef- 
fective Hamiltonian for the junction undergoing a phase 
slip and a quadratic Hamiltonian for the rest. 

The transformation of ri [Eq. (1)] to ri c n + "Hharm 
[Eq. (10) and Eq. (11)] to solve the first two of the three 
difficulties explained in the main text is most easily ob- 
tained via the Lagrangian C. For Cq — 0, a Legendre 
transformation yields 



n=0 



N-l 
n=0 



1 — COS 



S 

N 



(17) 



Note that we labeled the variable 9m by 9q for conve- 
nience. Because only N — 1 of the N phase difference 
variables 6i are independent [cf. Eq. (2)], the variable 
6>jv_i can be replaced by 



N-2 



9n-i = — 2, @n mod 2-7T . 



(18) 



n=0 



The phase differences enter Eq. (fT7j) only modulo 27T. 
Therefore, Eq. JTS]) leads to 



C 



K 2 C 
~8e Y 

h 2 C 



N-2 

E 

n=0 



9i - E.. 



'N-2 

E< 

\71=0 



JV-2 

E 



Ej 




We now focus on the case $ = $0/2, where S = tt 
and the lowest two PCS are \m = 0) and \m = —1). As 
mentioned in the main text, we represent the \m = 0) 
PCS by the minimum {9i — 0}. Let us first assume that 
the phase slip bringing the system into the \m = —1) 
PCS happens in junction no £ {0, . . . , N — 2}; the case 
no = N — 1 will be discussed at the end. We begin with 
a rotation of the coordinate system {0$, ... ,9^-2} -^ 
{9$, . . . , 9 N-2}, which relabels the junction no as junction 



e k = e k+no , fcG{l,...,JV-2-no}, 
Ok = 6 k +n a -N+i , k€{N-l-n ,...,N 



2}. 



(20) 

Next, a new set of variables {#o, • • • ,0n-2} -^ 
{yo,.-.,yN-2} is constructed by choosing yo along the 
phase slip direction and completing the change of vari- 
ables by a Gram-Schmidt process: 



1 



yo 



2/1 



2/n 



ViV 2 - AT - 1 



V^ 



VN 2 - N - 1 



(JV - l)0o 



iV- 1 



JV-2 

-E< 

j=i 

N-2 



VAT 



l * 



7V-2 ^ 
i=i 

1 



AT-2 



VTT 



y^T^nVN 



1 f^ 



(21) 

for n G {2, . . . , A^ — 2}. These variables have the desired 
property that the two current minima are related to each 
other by the shift of one coordinate (yo) only, while all 
other coordinate axes are orthogonal to it: 







o 



9„ = 2nS n . 



2tt 
~N 



o 



Vi : 

yo 

Vi 



o 



2tt\Jn 2 - N - \/N 

o »e{i,...,JV-2}' 

(22) 
After making a harmonic approximation of the poten- 
tial energies in Eq. (fT9|) for all 9 n ^ no (the cosine term 
must be kept for the variable 9 na = 9q), the resulting 
Lagrangian contains terms that couple different yi . This 
coupling is undesirable because it makes it difficult to 
evaluate the tunneling amplitudes. 

The above problem may be resolved by making an ad- 
ditional linear transformation. A coordinate system that 
satisfies the idea of Eq. (f22|) . while at the same time lead- 
ing to a completely decoupled Hamiltonian with the co- 
sine potential of junction no preserved, can be found to 
be 

\N-l)y + VN-2 yi ] :W„. 



Xl 



VN 2 - AT - 1 



VN 2 - N - 1 y/N -2 - 

x 2 = 7= = Vl = 



VN^T N ^\ 



VN 

2M-1 , 



n e {3, 



Vn- 

,N- 



1 
1} 



Vn- 



2 ^i 



After this transformation, Eq. (|22|) changes to 
9 n = ^ Xi = 



(23) 



9WS _^ 



^- 



xi = 2tt(N - l)/N 

Xi = ie{2,...,N-l} 



(24) 



&v( xi ) 




FIG. 4: Effective potential energy seen by the junction under- 
going a phase slip. Shown is the case N = 50, together with 
the harmonic approximations (dashed lines) about Xi — 
(PCS \m = 0)) and si = 2tt(N - 1)/N (PCS \m = -1)). 



and the Lagrangian reads 



2 AT- 



JV-l 



- if - V(ari) + -M J2 (il - <"X) , 

(25) 



n=2 



where (cf. Fig. @| 



Vfsi 



£.j 



1 — cos I x 



TV 



1 



■K 

A I ~ N Xl + 2(N - 1) 



(26) 

is the effective potential seen by the junction undergoing 
the phase slip. In Eq. (|25|) . M — h 2 /(8E c ) and uj p is 
defined in the main text [Eq. (6)]. Note that Eq. (|25|) 
contains no terms that couple different x,. From there, 
H —t H c ff + "Hharm follows immediately, with Eq. (10) 
and Eq. (11). 

For completeness, we provide the inverse transforma- 
tion of Eq. (f23| , which is necessary to arrive at the results 
presented in the main text: 



A) = xi 



1 



1 



■Xl 



E 



N - 1 VN~^2 VN~T 

1 



■r'2 



f-^VN^jVN-j-1 



VN - n - 2 

x i +1 + / at T Xn+2 

VN — n — 1 



(27) 

for n € {1, . . . , AT — 2}. In addition, we note the useful 
identity 



N-2 



'N-2 



n=0 V 3=0 / 



A r 



iV-1 



N ■ 



i*? + £ 



(28) 



We now discuss the case no = N — 1. The above 
changes of variables do not work for this case. Never- 
theless, following the same ideas still leads to a suitable 
set of new variables. It turns out that real and imaginary 



parts of the discrete Fourier transform of {0 n }, 



Xk+l 



with 



-</F=TRe(&) 

\/2Refe) 
Re(&) 

V^Im^fe-ArJ 



fc = 

fce{i,.. 

fc = iv ft 
k = N h 
k e {N h 






.,^-1} 
iV even 
AT odd 
M,...,iV-2} 

(29) 



(30) 



and 



& = 



1 



iV-2 

En Jw^i^ 

v n— 



fce{o,...,JV-2} ; 



(31) 

lead (upon replacing x n by £ n ) to exactly the same La- 
grangian [Eq. (|25[) ]. and with it to the same Hamiltonian 
H e « + Hharm, while respecting Eq. (|24]l. For these vari- 
ables the inverse transformation is essentially given by 



AT-2 



^^to 



Cfce" 



(32) 



and Eq. (|28jl still holds. Even though the contribution 
of junction TV — 1 might seem different at first, it leads 
to exactly the same effective Hamiltonian as the other 
junctions, as expected from symmetry arguments. 



Energy-splittings in presence of random charges 

In this section we derive the distribution of the energy 
splittings in the one-plasmon subspace, and discuss the 
existence of a "zero mode" in the limit e 3> |z/°'|. 

The energy-splittings are given by 2\fi n \ 1 ' 2 , where \i n 
(n £ {1,. . . ,N — 1}) are the eigenvalues of the hermitian 
matrix M = (z/ 1 ))'' V' 1 ) . Using Eq. (12), it is straightfor- 
ward to obtain Mu> in the eigenbasis of plasmon modes; 
however, the resulting expression is cumbersome and will 
not be shown here. The diagonalization of M becomes 
significantly simpler after Fourier transforming to (junc- 
tion) coordinate space: 



M n 



1 



JV-1 

at-tE' 
i,i>=i 



■N-i e 



'Mn 



(33) 



This results in M nn , = M^, + M^„ where 



M 



(») 



,(0)l 2 



e 2 +2eRe( I /' »z;) 



M^\ = ^Re 



.2 



^ Uj (l --;,-<•) -jj<Zn>, (34) 



and we have neglected subleading terms in 1/N. Since 
the matrix elements of M ' 1 ' are N times smaller than 
those of M(°\ we can at first approximate M nn > ~ M. 
This yields 



(0) 



f 



(o) 



/(0)|2 



e 2 + 2eRe(^z;), 



(35) 



which is the result quoted above Eq. (16) in the main 
text. After noticing that \v^°'\ and z n are approximately 
independent distributed for A^ ^ 1, we plug Eq. (f35|) 
in the standard definition of the conditional probabil- 
ity [1] and straightforwardly arrive at Eq. (16). The 
corresponding eigenfunctions are spatially localized, be- 
cause (i) M„ n / ex S n:rl ', and (ii) y/ n are typically non- 
degenerate (due to random charges). This implies that 
plasmon modes get spatially localized (at Co = 0) in 
presence of phase slips and random charges. 

So far we have neglected M 1 - 1 ' . Although its matrix el- 
ements are N times smaller than those of M^°\ there are 
N times as many of them. This suggests that the effect 
of M' 1 ) on /i„ can be relevant, at least for some of the 
eigenvalues. Here we will analyze this issue analytically 
in the regimes e 3> l^ - 1 ! and e <C |i/°)|. We will argue 
that Eq. (16) remains a good approximation for (at most) 
all-but-one eigenvalues. This statement is in agreement 
with what we find from a numerical diagonalization of 
(V(i))Vi) f or arbitrary |i/(°)|/e. 

Let us begin from the case e 3> |^°^|, which is relevant 
for some experiments (e.g. e ~ 10(|z^°)|) in Ref. [2|). In 
this regime, we may approximate 



M 



(i) 



N 



(36) 



(i) 



M. 



(0) 



Due to the separability of M, 

diagonalized analytically. The eigenvalue equation reads 



M™, can be 



E 



M (0) , + M 



(1) 



(A) 



which can be rewritten as 



i>. 



w 



t 



z n c\ 



MaV'.. 



(A) 



where A is the eigenvalue label and 



<-x 



Y,Zn>^. 



(37) 



(38) 



(39) 



Substituting Eq. {35} in Eq. (O and using \z n \ 2 = 1, we 
get 



e V 1 



M^ - fi x 



If c\ 7^ 0, this gives 



JV-l 



'■sE 



A^ ^ | r ,(o)|2 + e 2 + 2e Re(i/(o)2*)- / u A ' 



(40) 



(41) 
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The right hand side of this equation is a sum of N — 1 
terms with poles at p\ = p x . For a given configuration 
of z n and for N ^> 1, a graphical solution of Eq. (J4TJ) 



shows that N — 2 eigenvalues, given by /ia — p x (A € 
{2, . . . ,N — 1}), are densely packed between (e — |i/°) |) 2 
and (e + I*/ - 1 !) 2 . Their respective eigenmodes coincide 
with those of M™, i.e. they are spatially localized. In 
summary, Eq. (16) remains valid for N — 2 eigenvalues. 
However, the graphical representation also shows the 
existence of a single eigenvalue p\ such that p\ < (e — 
|;/°)|) 2 . The appearance of this eigenvalue outside the 
dense distribution formed by all the other eigenvalues 
is mathematically analogous to what happens in the 
Cooper problem of superconductivity [3| . There, the low- 
est eigenvalue (equal to the binding energy of a Cooper 
pair) lies below a quasicontinuum energy spectrum for 
pairs of unbound electrons. 



Anticipating that ^! <C |j4°'| 2 - 
Eq. (gTJ) as 



- e 2 — 2e I z4 ' I , we rewrite 



In the above derivation we have assumed c^ ^ 0. If 
c\ = 0, we are constrained to a N — 2 -dimensional plane 
(orthogonal to the vector corresponding to A = 1) and as 
expected we get p\ = p x . 

Finally, let us briefly consider the opposite case 
|i/°)| ^> e. In this regime, we may approximate 



M. 



(i) 



2e 



^Re[^)(l 



N 



<>)] 



(45) 



and proceed to analyze the corresponding eigenvalue 
equation much like for the e ^$> |i/°)| case. This eigen- 
value equation is less simple than the one above; it reads 



- 1 = 



I-i dx -c 



p( x ) 



rl j p(x)x(l—x) 

J-i ax x-i-ix 



1 + f 1 dx p{x){1 ~ x) 
1 + J-i ux x-l-£\ 



j\dx 



p(x)x 



X-1-& 



(46) 



where fr = (/* A - e 2 - |^°)| 2 - 2e|i/°)|)/(2e|z/°)|), and 



'4E 



| I/ (0)|2 +e 2_ Ml 



2eRe(i/°>z*) 

I/ (0)|2 +e 2_ /il 



A simple calculation then shows that 



/'■i 



,(°)| 2 



2 £ |t/°)| 2 

A^ ! 



(42) 



(43) 



where we have used z/ ' = v J2 n z n . Thus, 2|^i| 1 / 2 
agrees approximately with the ground-state energy- 
splitting. Since e 3> |iA°)|, /ii <C Ma>i and thus the 
A = 1 eigenvalue may be loosely referred to as a "zero 
mode" . Its eigenfunction is 



v 



(A=l) 



oc 



M (0) 

lvinn 



(44) 



Mi 



modulo a n— independent normalization factor. For e 3> 
|^(°^|, the weight \tpn | 2 is approximately independent of 
n and hence the "zero mode" is completely delocalized 
onto all the junctions in the ring. As e/|i/°)| becomes 
smaller, \p\^\ — Pi\ decreases, localization effects start 
to become visible in ?/4 , and the term "zero mode" grad- 
ually becomes meaningless. 



1 



JV-l 



<°( x ) = A? J2 S(x-Re(z n )). 



(47) 



The function p(x) is random. In the derivation of 
Eq. (|46)) . we have made a phase rotation so that z4°) -4 
|iA°)|. We can gain some insight into the solution of 
Eq. (J46]) by replacing p(x) -> (p{x)) ~ (l/7r)(l-a; 2 ) _1 / 2 , 
upon which Eq. (|46[) can be solved numerically. The so- 
lution shows that p\ ~ /i^ is a good approximation for 
all N — 1 eigenvalues (provided that N 3> 1). In other 
words, there is no eigenvalue that lies outside the distri- 
bution of Eq. (16). Thus, the appearance of a delocalized 
"zero mode" is exclusive to the regime e ^> \u^°'\. 
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